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ABSTRACT 

Sagittarius A* harbors the supermassive black hole that lies at the dynamical center 
of our Galaxy. Sagittarius A* spends most of its time in a low luminosity emission state 
but flares frequently in the infrared and X-ray, increasing up to a few hundred fold in 
brightness for up to a few hours at a time. The physical processes giving rise to the 
X-ray flares are uncertain. Here we report the detection with the NuSTAR observatory 
in Summer and Fall 2012 of four low to medium amplitude X-ray flares to energies up 
to 79 keV. For the first time, we clearly see that the power-law spectrum of Sagittarius 
A* X-ray flares extends to high energy, with no evidence for a cut off. Although the 
photon index of the absorbed power-law fits are in agreement with past observations, 
we find a difference between the photon index of two of the flares (significant at the 95% 
confidence level). The spectra of the two brightest flares (~55 times quiescence in the 2- 
10 keV band) are compared to simple physical models in an attempt to identify the main 
X-ray emission mechanism, but the data do not allow us to significantly discriminate 
between them. However, we confirm the previous finding that the parameters obtained 
with synchrotron models are, for the X-ray emission, physically more reasonable than 
those obtained with inverse-Compton models. One flare exhibits large and rapid (< 

100 s) variability, which, considering the total energy radiated, constrains the location 
of the flaring region to be within ~10 Schwarzschild radii of the black hole. 
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1. Introduction 

Sagittarius A* (Sgr A* ), at a distance of ~ 8 kpc, is a ~ 4 x 10 6 solar mass supermassive black 
hole (SMBH) that marks the dynamical center of the Milky Way (Ghez et al. 2008; Gillessen et al. 

2009). Although Sgr A* is a fairly bright radio and submillimeter (sub-mrn) source, it is not detected 
in the optical and UV due to over 30 magnitudes of visual extinction, and it is very dim in X-rays. 

Its bolometric luminosity L^ 0 \ is only 550 times that of the Sun, corresponding to ~ 10~ 9 times 
the Eddington luminosity (Narayan et al. 1998). Most of this is radiated in the sub-mrn, with 
the X-ray luminosity being 3.6 x 10 33 erg s" 1 (Nowak et al. 2012) in the 2-10 keV energy band 
(~ 2x 10~ 3 Abol)- This low luminosity implies that matter falls onto the black hole in a radiatively- 
inefficient manner, perhaps in a hot, geometrically thick, optically thin flow (Narayan et al. 1998; 

Yuan, Quataert & Narayan 2003; Wang et al. 2013). 

A few times per day in the near infrared (NIR) and about once per day in X-rays, Sgr A* flares 
up to a few hundred times its quiescent level for intervals lasting up to a few hours, with the 
strong flares being less common than the weak ones. The NIR emission is strongly polarized 
(Eckart et al. 2006b), indicating a synchrotron origin, and therefore a population of relativistic 
electrons. The source of the X-ray emission is uncertain; possibilities include both synchrotron ra- 
diation with a cooling break (SB, Yuan, Quataert &; Narayan 2003, 2004; Dodds-Eden et al. 2009) 
and inverse Compton (IC) up-scattering of lower energy photons by energetic electrons. Two IC 
scenarios involve the NIR emitting electrons up-scattering either ambient sub-mm photons (exter- 
nal Compton, EC, Markoff et al. 2001; Eckart et al. 2004; Yusef-Zadeh et al. 2006a) or the NIR 
synchrotron emission itself (synchrotron self-Compton, SSC, e.g. Markoff et al. 2001; Eckart et al. 

2006a; Marrone et al. 2008). A third IC scenario involves NIR flare photons up-scattered by the 
electrons radiating in the sub-mm domain (Yusef-Zadeh et al. 2009, 2012). 

The physical cause of the episodic particle acceleration is not understood, and suggestions 
include a hot spot in the accretion flow (due, for example, to a magnetic reconnection event, see 
e.g. Dodds-Eden et al. 2010), enhanced mass accretion (Tagger & Melia 2006), magnetohydrody- 
namic turbulence or shocks in the inner accretion flow (Yuan, Quataert Sz Narayan 2003), a jet 
(Markoff et al. 2001), or episodic outflow triggered by magnetic reconnection (Yuan et al. 2009). 

Tidal disruption of asteroids has also been proposed as possible origin of the flares (Cadez, Calvani &; Kostic 
2008; Kostic et al. 2009; Zubovas, Nayakshin & Markoff 2012). 

Numerous flares have now been observed in X-rays (Baganoff et al. 2001; Porquet et al. 2003, 

2008; Degenaar et al. 2013; Neilsen et al. 2013), in NIR (Genzel et al. 2003; Witzel et al. 2012), 
and in both bands simultaneously (Eckart et al. 2004; Yusef-Zadeh et al. 2006a; Eckart et al. 2006a; 
Marrone et al. 2008; Dodds-Eden et al. 2009; Trap et al. 2011; Eckart et al. 2012). Above 10 keV, 
only upper limits exist due to a combination of limited instrumental sensitivity and spatial reso- 
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lution (Yusef-Zadeh et al. 2009; Trap et al. 2011). Short, high-amplitude temporal substructures, 
with variability timescales as short as 47 s, have been seen during NIR flares (Dodds-Eden et al. 
2009), indicating a compact origin for this emission component. In the X-ray, variability on 
100-200 s timescales has been reported, although at low relative amplitude (Porquet et al. 2003; 
Nowak et al. 2012). The fastest timescale observed to-date for significant (more than a factor of a 
few) X-ray variability is 600 s (Baganoff et al. 2001). 

We report here on NuSTAR observations of four flares detected in Summer and Fall 2012. This 
paper is organized as follows. Section 2 details the observations, and section 3 presents the method 
to search for flaring activity. In section 4, we quantify the variability amplitude and time scale 
exhibited during one of the flares. Section 5 describes the spectral analysis of the flares. Then, in 
section 6, we introduce the models that are compared to the SEDs of the two brightest flares, and 
the results are discussed in section 7. Finally, in section 8, we derive a constraint on the location 
of the flaring region. 


2. Observations 

The NuSTAR high-energy X-ray observatory (Harrison et al. 2013) observed the Galactic cen- 
ter three times in Summer and Fall 2012 as part of a coordinated campaign with the Chandra 
and Keck observatories: from July 20, 2:11 to July 23, 19:11, from August 04, 07:56 to August 
06, 01:06, and from October 16, 18:31 to October 18, 01:26 (UT time). Flaring activity from Sgr 
A* was detected by NuSTAR on July 20 and 21, and on October 17. Four flares were detected, 
which are referred to hereafter by J20, J21_l, J21_2 and 017. 

Unfortunately, only 017 was simultaneously observed by Chandra , and none of the flares were 
covered by Keck. We report in this paper on the NuSTAR data only. The Chandra data for the 
017 flare will be presented in a future work. 


3. Flare search: Bayesian block analysis 

The data analysis was done with the NuSTAR Data Analysis Software ( NuSTARDAS ) vl.1.1. 
We first screened the data for South Atlantic Anomaly (SAA) passages and other flaring particle 
backgrounds. Then, we used an extraction region with a 50 7/ radius (encompassing ~ 68% of the 
PSF) centered on the radio position of Sgr A* (Reid &; Brunthaler 2004) to extract the events from 
focal plane modules (FPM) A and B in the full energy band (3 - 79 keV), and we performed a 
Bayesian block analysis of the combined light curves to search for flares. This analysis is applied 
on the raw events of both FPMs after the data gaps due to Earth occultation and SAA passages 
are suppressed. We correct the exposure by applying a coefficient to each event corresponding 
to the deadtime at the time of the event. The Bayesian block analysis assumes that the light 
curve can be modeled by a succession of constant rate blocks. The rate within each block is 
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Fig. 1. — The Sgr A* field (4.5' x 4.5') during flare J21_2 (left panel) and before the flare (right 
panel) as seen by NuSTAR’s focal plane module A. The images are made with 3 to 79 keV photons 
without background subtraction and an exposure time of 3000 s for each panel. It is smoothed with 
a 3-pixel radius Gaussian function. The image is oriented in equatorial coordinates with North to 
the top and East to the left, with a logarithmic color scale. The grid shows Galactic coordinates, the 
cross shows the Reid et al. (2004) radio position of Sgr A* , and the small circle (dash line) shows 
the 3 -o NuSTAR error circle. The large circle shows our 50" radius extraction region, centered on 
the radio position. 

simply calculated based on the number of events it contains and its duration. The algorithm to 
find the optimal location of the change points is described in Scargle et al. (2013). The number 
of change points is affected by two input parameters: the false positive rate fpr and the prior 
estimate of the number of change points, n cp _ pr i or . We used fpr = 0.01 and a geometric prior: 
^cp.prior = 4 — log(/pr/(0.0136 IV 0 - 478 )), where N is the total number of events (Scargle et al. 
2013). 

We investigated the actual false positive rate by performing simulations with synthetic datasets 
having a constant input rate (equal to the average rate of our dataset) and following a Poisson 
distribution for the event arrival time intervals. We obtain identical results for both datasets (the 
July and October light curves have 52,109 and 28,067 events, respectively). We find that spurious 
peaks (a peak is made of two change points) are detected 0.4% of the time (1000 trials performed). 
See Appendix A for more details on the sensitivity of the search. 

This analysis led to the detection of four flares as listed in Table 1. The flare blocks are 
used to define the flaring times. The quiescent time is defined independently for each observation 
and includes the time remaining after removing the flares and the flaring background events. The 
quiescent time for the July and October observations amounts to 53.5 ks and 41.8 ks, respectively. 
Figure 2 shows the July and October portion of the light curve where the flares were detected, and 
Figure 3 shows close-up views of the flare light curves. NuSTAR is in low-Earth orbit, and the gaps 
in the data are due to Earth occupation and passages through the SAA. 
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The position of the flaring source is aj 2 ooo = 17 fe 45 m 40.4 s , <5j2ooo = — 29°00'31.04' / , with an 
error circle of 10. 5" radius (3-<r confidence level), including a systematic component estimated to be 
at the 5" level (Figure 1). The position is consistent with the radio location of Sgr A* (5. 6" offset; 
Reid <fc Brunthaler 2004). The association with Sgr A* is clear; no variable source was detected 
by Chandra within the error circle during the July 20-23 and October 17 observing campaigns (F. 
Baganoff, private communication), and the characteristics of the flares are similar to those seen 
in prior low-energy X-ray observations. No other source is known to produce X-ray flares with 
luminosity of the order of 10 3 ° erg/s, a duration of the order of an hour, and a power-law spectrum 
with photon index around 2. 


4. July flare variability 

Substructures are clearly visible in the light curves of flares J21_2 and 017. In particular, on 
2012 July 21, there is a rise at 6:21 UT and a decay at 6:51 UT that both happen within one 
100-s bin (Figure 3). These features are visible in both FPM A and FPM B (see Appendix B). 
To estimate the amplitude of these variations, it is necessary to understand that the light curve is 
offset due to the unresolved high-energy emission that is present within our 50 w radius extraction 
region, and that this diffuse emission dominates the count rate at all times except during flares 
J21_2 and 017. This baseline count rate is measured to be 0.585 ± 0.011 cts/s, as averaged using 
the blocks preceding and following flare J21_2. 

We first focus of the rise at 6:21 UT. Estimating its amplitude is not straightforward as the 
bottom point is at the level of the quiescent flux (a statistical fluctuation makes it appear below the 
quiescent level, visible in the bottom left panel of Fig. 3), consistent with a non-detection during 
the dip that separates sub-flares J21_2a and J21_2b. Near the time of the rise, the bottom and top 
bins show rates of 0.55lQQg cts/s and 0.92 +q'Jq cts/s, respectively. We use the upper 1 -a error bar 
as the upper limit of the flare count rate during the minimum, which leads to a lower limit on the 
rise factor of 3.8 ± 1.1. During the 6:51 UT decay, the count rate drops from 1.22+Q12 cts/s to 
0.771 q'q 9 cts/s. Subtracting the baseline count rate from both these values, we find a decay factor 
of 3.4 ± 1.8. 


5. Spectral analysis 

The flare and quiescent spectra are extracted from the same 5CU radius region centered on 
the radio position of Sgr A* . The analysis is done with the Interactive Spectral Interpretation 
System (ISIS, Houck h Denicola 2000a) by setting the non- flare spectra as background to the flare 
spectra. The model TBabs x PEGpwrlw (i.e., an absorbed power-law) is used to fit the data, with 
Verner et al. (1996) atomic cross sections and Wilms, Allen &; McCray (2000) abundances for the 
interstellar absorption. 
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Fig. 2. — Light curves from a 50" radius region centered on the radio position of Sgr A* in the 
3-79 keV energy band for NuSTAK s modules A and B combined. The top panel shows a portion 
of the July observation, and the bottom one shows a portion of the October observation. In both 
panels, the grey histogram has 100 s time resolution and shows 1-a upper and lower count rate limit. 
The red lines show the Bayesian block representation of this light curve (Scargle et al. 2013), with 
1-a rate limits as dotted lines. Instrumental background, estimated to contribute < 5 x 10 -3 cts/s 
in this energy band for this extraction region, has not been subtracted in this plot. The diffuse 
emission from the Sgr A* complex is the main component to the ~ 0.6 cts/s that constitutes the 
baseline of the light curve. The flares IDs are shown in the plot. 


Table 1. Flares IDs, start and stop times (UT), duration, and significance (for FPM A and B 

combined) . 


Flare ID 

UT Date 

Start 

Stop 

Duration (s) 

Level (a) 

J20 

2012 July 20 

12:15:21 a 

12:30:39 

920 a 

4.8 

J21.1 

2012 July 21 

01:45:15 b 

02:06:14 b 

1238 b 

7.5 

J21.2 

2012 July 21 

06:01:12 a 

06:52:50 

3099 a 

19.7 

017 

2012 October 17 

19:50:08 b 

20:10:56 b 

1249 b 

20.2 


a This flare was truncated at the beginning by Earth occultation. 

b This flare was truncated at the beginning by a passage in the SAA and at the 
end by Earth occultation. 



Fig. 3. — Close up view of the light curves and associated Bayesian block representations for the 
four flares reported in this paper. Binning and error bars are identical to those in Figure 2. The 
shaded areas show the baseline count rate as defined by the bayesian blocks preceding and following 
each flare. 
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Fig. 4. — Counts spectrum of the four flares extracted from a SO^ radius circular region centered 
on the Sgr A* radio position: J20 (black crosses), J21_l (red triangles), J21_2 (green squares), and 
017 (blue plus signs). Each spectrum is the combination of data from modules A and B. The data 
points are shown with 1-cr error bars, and the solid lines show the best fit model. The lower panel 
shows the deviation from the model in units of standard deviation. 



Fig. 5. — Photon index (T) of the flares (left panel), and for a breakdown of flare J21_2 (right panel), 
with the letters r, p, and d referring to rise, peak and decay, respectively. In the left panel, the four 
flares are fit jointly by an absorbed power-law model with the column density tied (1.66jjo g? x 10 2 ^ 
cm -2 ). The breakdown of flare J21_2 uses a column density fixed to the value determined using 
the four flares. The errors bars show the 90% confidence range. 
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Some previous studies of the spectra of Sgr A* flares included a dust scattering correction 
Baganoff et al. (2001); Porquet et al. (2003, 2008); Nowak et al. (2012). The effect of dust is 
twofold; it contributes to absorption (through the photoelectric effect), which is accounted for 
by TBabs in our model, and it scatters X-rays, primarily in the forward direction, creating a halo. 
The scattering has two consequences: it removes photons from the direct line of sight (LOS) to the 
object, and it sends photons initially emitted out of the LOS towards the observer, creating the 
halo (Predehl et al. 2000). With an extraction region of 50 ,/ radius, we can assume that these two 
effects compensate each other, so we do not include the loss due to dust scattering in our model. 
We note that we find a column density and a photon index consistent with past studies, a cross 
check on this assumption. 

Figures 4 shows the spectra measured for the four flares with their best fit absorbed power-law 
model. The spectral binning is identical for each spectrum: spectra of FPM A and B are combined, 
and the spectral bins are grouped to reach at least 50 source + background counts and 3 a signal- 
to-noise significance in each group. The high energy bound of the last group is set to maximize the 
significance of the group. Finally, the last group is merged with the second to last if it does not 
reach 2 a significance. The last bin significance is 2.51 (8.52-21.00 keV), 2.79 (9.88-26.00 keV), 2.88 
(22.76-60.00 keV) and 3.77 (20.16-79.00 keV) for flares J20, J21_l, J21_2, and 017, respectively. 

We assume that the column density is constant from flare to flare, and we perform a joint fit of 
the four flares with the column density tied together. We find a column density of 1.66lto gi x 10 23 
cm -2 , with a reduced y 2 of 0.90 (52.21/58). For comparison, column densities measured during past 
flares range from 1.25 x 10 23 cm' 2 to 1.85 x 10 23 cm' 2 (converting published values to be consistent 
with Verner et al. cross sections and Wilms et al. abundances, Goldwurm et al. 2003; Porquet et al. 
2003; Belanger et al. 2005; Porquet et al. 2008; Nowak et al. 2012), being all consistent with each 
other within the uncertainties (which is consistent with our assumption of the column density 
being constant from flare to flare). Table 2 gives the best fit photon indices and associated fluxes 
and luminosities for each flare. The photon indices we find are consistent with those reported in 
previous work (Porquet et al. 2003; Belanger et al. 2005; Porquet et al. 2008; Degenaar et al. 2013; 
Nowak et al. 2012). 

Flare 017 residuals show possible evidence for a cut off around 20 keV. For a direct comparison 
between the spectra of flares J21_2 and 017, which have similar peak luminosity, we produced 
another 017 flare spectrum with an upper energy limit of 60 keV, and the highest energy point 
still falls below the power-law. However, fixing the column density to the value found for the four 
flares, an F-test shows that a cut off is not required by the data (probability of 90.2%). 

We split flare J21_2 into four phases to look for spectral evolution (see Figure 3). The spectral 
binning is made in the same way as described above, and the column density is fixed to 1.66 x 10 23 
cm -2 . Figure 5 shows the evolution of the spectral index during flare J21_2 (right panel), and the 
best fit photon index of the power-law for all four flares (left panel). There is no spectral evolution 
during flare J21_2, as T changes by significantly less than its error across the flare, which confirms 
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past results (Porquet et al. 2008; Nowak et al. 2012). 

However, the 90% confidence ranges of flares J21_l and 017 do not overlap. We investigate 
the significance of this difference by plotting the confidence contours of both photon indices (Figure 
6). The calculation is done similarly to the spectral analysis presented above: the four flares are 
fitted simultaneously with tied column density. We find that the two photon indices are different 
at the 95% level. We know, thanks to Chandra, that we caught flare 017 right at its peak (J. 
Neilsen, private communication), but, as previously noted, there should not be any bias due to 
spectral evolution. Although this result is only marginally significant, we note that it is the first 
time such a difference has been reported between two individual flares. Degenaar et al. (2013) 
noticed a low-significance difference by comparing one bright flare to the average of five weak flares. 
Interestingly, the bright flare was softer (T 2 -iokeV,peak = (3.1 ± 0.5) x 10 35 erg s -1 , T = 3.0 ± 0.8) 
than the average of the five weaker ones, which is opposite to what we observe. This discrepancy 
can be interpreted as a piece of evidence that Sgr A* X-ray flares can have a range of photon 
indices without being correlated to luminosity, keeping in mind that both the present result and 
that of Degenaar et al. (2013) have low significance. We note that magnetic reconnection is known 
to create MHD turbulence that produces stochastic acceleration, which can result in a range of 
particle distribution slopes (Zharkova et al. 2011). This could explain flares of various photon 
indices. 


6. SEDs and one-zone static models 

An important step in understanding the production of the flares is to determine the underlying 
X-ray emission mechanism. Synchrotron emission produces a spectrum following a power-law. In 
the SB case, where the X-ray flare is produced by the same electron population that emits the NIR 
flare, a spectral break occurs in the UV or optical due to fast electron cooling (Dodds-Eden et al. 
2009). Alternatively, we consider two IC scenarios, SSC and EC, both involving the electrons 
producing the NIR synchrotron emission flare (the third IC scenario mentioned in the introduc- 
tion is not implemented here as it produces a similar X-ray spectrum to the other two). All 
three mechanisms could originate either in a jet or in an accretion flow. It is not our ambition to 
perform detailed modeling of the flaring emission. Instead, we want to identify the dominant mech- 
anism by using single-component, one-zone (all parameters are constant across the emission region), 
and static models following Dodds-Eden et al. (2009). In previous work, Eckart et al. (2012) and 
Yusef-Zadeh et al. (2009) used more sophisticated dynamical models to fit multi-wavelength light 
curves. Our simpler approach is appropriate given our limited bandpass (i.e. X-rays only). The 
constraints that we derive are reasonably robust, although some of the inferred parameters can 
change in more sophisticated flare scenarios (e.g., the magnetic field strength, Eckart et al. 2012). 

We calculate the X-ray spectra of the IC models according to the equations 18 to 21 of 
Chiaberge & Ghisellini (1999) (see also Rybicki & Lightman 1979). The free parameters for the 
SSC scenario are the electron density n e , the size of the emitting region R, the magnetic field 
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Fig. 6. — Confidence contours for the photon indices of flares J21_l and 017, determined while 
fitting simultaneously the four flares with tied column density. The four contours show the 68.3%, 
90%, 95% and 99% confidence levels (solid line, dash line, dot line, and dash-dot line, respectively). 
The straight line shows the values where the photon indices of both flares are equal. One can see 
that the photon index of flares 017 and J21_l are different at the 95% level. 


Table 2. Photon index, fluxes and luminosities for each flare in the 2-10 keV and 3-79 keV 

energy ranges. 


Flare ID 

r 

ZTiabs rpunabs 

* 2-10 * 2—10 

( 10 -11 ergcm - 2 s -1 ) 

t unabs 
*^2-10 

( 10 35 ergs -1 ) 

77 abs runabs 

r 3-79 *3-79 

( 10 -11 ergcm - 2 s -1 ) 

t unabs 
*^3-79 

( 10 35 ergs -1 ) 

Strength 

J.20 

2 69 +0 ' 98 
z,.u»_ 0 84 

0-38±g;iJ 

0 «<5+ 0 - 62 
U.SO-o 39 

0 f)>5+ 0 - 48 
uoo -0.30 

0-TOtg;gg 

0 95 + 0 - 50 

U.yo_o 34 

n 7q+0.38 
u - ‘°-0.26 

18.i! 13 4 3 

J21.1 

2 84" 1 ” 0 ' 64 
Z>O ^-0.54 

o-5i±8:ii 

1 20”*”°' 63 
±,zu -0.41 

0 92”*”°' 49 
u ’ yz -0.31 

0-86tg;i* 

1 20"*”°' 33 
lzu -0.30 

0 92”*”°' 25 
u,y -0.23 

25.5l 13 9 7 

J21.2 

9 qq+0.24 

z.zO—o .22 

n 09 + 0.09 

U.84_o 09 

1 co+0.43 
i.oo _ 0 32 

1 oc+0.33 
t.zO-024 

9 99+O.37 

z,zz -0.33 

9 79+O.34 
z - ' z -0.32 

9 ns+0-26 

Z.Uo_o.24 

34.7±|;§ 

017 

9 04"*"°' 22 
z - U4 -0.20 

1 qq+014 
±.OO _ 0 13 

9 1-9+O.64 

l-93±g;$? 

4 41 + 0 - 72 

4 - 41 -0.65 

c 1 O+0.65 
o.io _ 0-61 

q Q7+0.50 
1 -0.47 

53.6±ig;| 

J21_2a 

1 90 + °' 48 

i.yu_o 4i 

0-48±C41 

p 

bo 

00 

1 + 
0 0 
to CO 
00 ^ 

0 . 68 l °;=6 

1 RrTO. 84 
l-oo_0.65 

9 1 ec+O- 8 ! 
z * * <J — 0.63 

1.0°_o 48 

18.71 7 ; 3 

J21_2b 

9 qi +0.19 
Z, °- L — 0.18 

1 02" 1 ” 0 ' 11 
± * uz _0.11 

2 06”*”°' 28 
Z.UO_Q 27 

1 t;O + 0.22 

l.OO-o 20 

2 5R - * -0 - 40 
z ' oo -0.36 

q 1 q+0.39 
*°— 0.36 

2 43+ 030 
* ,40 _Q,27 

43.8l 7 ;t 

J21_2b peak 

2 44 - *" 0,33 
z,44 -0.29 

1 97+O.2I 
*- z * -0.20 

2 66”*"°' 60 
ZOO -0.54 

2 f)3 + 0 - 46 
Z.UO_ 0 41 

9 qc+0.72 

z.oo_q go 

3 fi 4 + 0 - 70 

9.04-0 62 

9 7Q+0.53 

z - ‘ y -0.47 

56.6±i|;S 


Note. — The second column gives the best fit photon index T. The fluxes are determined using the cflux convolution 
model. The column density of 1.66 ^q g? x 10 23 cm -2 was determined by jointly fitting flares J20, J21_l, J21_2 and 017 with 
the column density tied together. The column density for the breakdown of flare J21_2 (in the lower part of the table) was 
fixed to 1.66 x 10 23 cm -2 . Both absorbed (noted abs) and unabsorbed fluxes (noted unabs) are reported. The luminosity 
calculations assume a distance of 8 kpc and isotropic emission. The strength (rightmost column) is defined as the ratio of 
the 2-10 keV unabsorbed flux to the 2-10 keV unabsorbed quiescent flux, 0.47^q'q3 x 10 — 12 erg cm _2 s _1 , as reported in 
Nowak et al. (2012). All errors are at the 90% confidence level. 
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strength B , and the temperature of the electrons T e . Assuming a thermal distribution of electrons, 
the spectrum due to synchrotron emission is calculated, which is then up-scattered by the same 
population of electrons. Double up-scaterring is not accounted for in the model. 


The EC scenario has additional free parameters, but we only vary those that affect the X-ray 
emission: n e , T e , R, and the size of the sub-rnm emitting region, -R S ub-mm- The sub-mm spectrum 
is taken from the Yuan, Quataert & Narayan (2003) model fit to the quiescent Sgr A* spectrum. 
Again, only one up-scattering is accounted for, and due to the absence of NIR constraint, no 
synchrotron emission is calculated (its intensity would only depend on B). 


The free parameters for the SB case are the index of the power-law distribution of the electron 
energy spectrum p, the magnetic field B, the total number of electrons N e (which strongly depends 
on the minimum Lorentz factor 7 m in)> and the injection timescale Ti n j (time for N e electrons to be 
accelerated). The cooling break creates the following electron distribution: 


71(7) 


7 p 7min < 7 < 7 c 
7 -(p+i) 7 C < 7 < 7 max 


where 7 C is the Lorentz factor of the cooling break. An electron population with power-law index 
p creates a power-law radiation spectrum with index (p — l)/2. 7 C is given by the balance between 
injection of energetic electrons and synchrotron cooling. We estimate r; n j based on the variability 
timescale observed during flare J21_2, 100 s. The cooling time, r coo i, depends on the magnetic 
field strength and on the energy of the electrons (the shortest cooling time is given by 7 max ). 7c is 
determined by the condition (Begelman, Blandford & Rees 1984): 


Tcool — Tnj 


3 8Trm e c 

4 ctt 7c B 2 


100 , 


where m e is the electron mass, c is the speed of light, or is the Thomson scattering cross-section for 
electrons, and j3 e is the ratio of the magnetic energy density to the thermal energy density. Solving 
this equation gives y c = 8 x 10 6 x B~ 2 (assuming equipartition, j3 e = 1). This induces a break in 
the synchrotron spectrum at the frequency 




0.457^ v gyT sin(a) ~ 6.3 x 10 7, -y 2 



6.4 x 10 14 



Hz 


where v gyr is the electron gyration frequency and a is the angle between the electron motion and 
the magnetic field direction (assumed to be 90° here). For flare J21_2, with an injection timescale 
of 100 s and an assumed value of 50 G for the magnetic field strength, the break is in the visible. 


The three aforementioned models were fit to the data using the chi-squared minimization 
method. Figure 7 compares the spectral energy distributions (SEDs) of flares J21_2 and 017 to 
the three models. The SEDs were calculated by determining fluxes and errors for modules A and 
B separately and then averaging them. A coarse binning was chosen for the sake of clarity (see 
Appendix C for the significance of the bins). We corrected the points for interstellar absorption as 
described in Appendix C. The best fit parameters for the three models are shown in Table 3. 
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Fig. 7. — Comparison of the NuSTAR SED for flare J21_2 (a,b) and 017 (c,d) to three one-zone 
models: SB (solid line), SSC (dashed line), and EC (dotted line). See text for details. The vertical 
lines of the data points show the 1 -a error bars. The lower panels show the residuals of the SB 
model. 


7. Implications of the SED measurements 

Work to date combining NIR and X-ray data does not definitively rule out any of the above 
mechanisms. The narrow bandpasses of Chandra and XMM-Newton and the strong absorption 
in the soft X-ray band have made it impossible to distinguish between a power-law spectrum, 
naturally explained by a synchrotron process, and a curved spectrum that would be characteristic 
of straightforward EC or SSC emission models. NuSTAR , with its larger bandpass, clearly shows 
that the spectrum extends to high-energy X-rays, but we are limited by statistics in the highest 
part of the energy range. The fit to flare J21_2 is somewhat conclusive with a rejection probability 
for the SB model of 23% versus 60% and 59% for the EC and SSC models, respectively. However, 
the fit to flare 017 is clearly inconclusive with no model standing out. 

We nonetheless note that the best fit parameters for the EC and SSC models are unrealistic, 
while the SB model gives reasonable ranges of electron density and magnetic field strength. The 
EC process requires an electron density several orders of magnitude higher than expected for the 
accretion flow around Sgr A* , and the high density required for the sub-mm photons leads to a very 
small sub-mm volume, over an order of magnitude smaller than the quiescent emission region as 
observed by VLBI (Fish et al. 2011). The SSC model needs even higher electron density along with 
extremely high magnetic field strengths. Thus, we reach identical conclusions to Dodds-Eden et al. 
(2009): we favor the SB model primarily based on physical arguments. 

Within the synchrotron picture, NuSTAR ' s high-energy detection implies either a higher mag- 
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Table 3. Best fit parameters for the EC, SSC and SB models of flares J21_2 and 017. 


Parameters 

EC J21_2 

EC 017 

SSC J21_2 

SSC 017 

SB J21_2 

SB 017 


n e (cm 3 ) 

4.2 x 10 7 

3.7 x 10 8 

3 x 10 10 

R (cm) 

4.6 x 10 10 

3.2 x 10 10 

5.4 x 10 9 

T^sub — mm (cm) 

1.1 x 10 11 

6.4 x 10 10 


B( G) 



3200 

Te (K) 

1.0 x 10 12 

1.4 x 10 12 

1.0 x 10 11 

Tinj(s) 




NIR flux (mJy) 



5.5 

X 2 / dof 

7.3/7 

5.4/7 

7.2/7 

CDF 

0.60 

0.39 

0.59 


3.4 xlO 10 > 2.9 — 1.9 x 10 3 >5.7-5.6xl0 2 


5.4 x 10 9 

< 3 x 10 12 

< 3 x 10 12 

2300 

2.7-24 

6.3-56 

1.4 x 10 11 

100 

100 

5.1 

1 - 27 

1-27 

4.8/7 

4.1/7 

4.4/7 

0.31 

0.23 

0.27 


Note. — The NIR flux is the value for the K s band (2.2 /im); in the case of the SSC model, the NIR flux 
is the prediction from the model. In the case of the SB model, it is simply the range of past observations 
(Trap et al. 2011). The EC model does not allow a NIR flux prediction, as it only depends on the magnetic 
field strength given the number of electrons and their temperature. In the SB case, the number of electrons is 
calculated for 7 min = 1000 and 7 max = 3 x 10 5 . A range of electron densities and magnetic field strengths are 
quoted, which produce the same X-ray flux but cover the range of NIR flux measured in past flares. The size 
of the emitting region does not play a role. The flux only depends on the total number of electrons. However, 
for consistency, we use the upper limit on the flaring region size that we derived for flare J21_2 to estimate 
the lower limit on the electron density (see Section 8). The last row, CDF, gives the cumulative distribution 
function for the y 2 distribution, given the x 2 value and the degree of freedom (dof). 
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netic field strength or a particle distribution with a higher maximum Lorentz factor (E max oc Bj^ &x ) 
than previously assumed. With a magnetic field strength of 50 G, flare J21_2 requires 7 max ~3 x 10 5 , 
which implies a synchrotron cooling time as short as ~ls for the most energetic X-ray emitting 
particles. Thus, continuous acceleration of high-energy particles is required to produce a ~3,000s 
flare, even if the magnetic field were more than an order of magnitude lower than the value given 
above. 

The rapid X-ray variability detected by NuSTAR also puts important constraints on flare 
models. In the SB model with a cooling break between the NIR and X-ray, the NIR variability 
is primarily due to changes in the magnetic field strength, while any X-ray variability is driven 
by sporadic injection of high-energy particles into the emission region (Dodds-Eden et al. 2010). 
Thus, the amplitude changes (rise and decay) within 100 s measured during flare J21_2 require that 
the particle injection rate can either drop or rise dramatically on a timescale of 100 s (see section 
4). Given that we observe these fast variations at the beginning and at the end of sub-flare J21_2b 
(separated by 30 minutes), we assume that the emitting region keeps the same size throughout the 
flare. For this given flare, this assumption puts constraints on the adiabatically expanding blob 
model for the X-ray emission (van der Laan 1966; Eckart et al. 2006a, 2012; Yusef-Zadeh et al. 
2006b, 2008, 2009). 

8. Location of the emitting region: energy balance calculation 

Regardless of the emission mechanism and of the scenario considered (disk or jet), we estimate 
the amount of available energy for producing flares using the magnetic energy density (B 8 * 2 * /8tt) and 
the volume of the flare’s emission region (R). Considering theoretical models that predict magnetic 
field strengths in the accretion flow and at the base of the jet, we compare the energy available at 
different distances from the SMBH to the fluence of the flare. 

The natural speed associated with the magnetic energy density is the Alfven speed, but in the 
inner radii of the SMBH, this speed is likely to be close to the speed of light c. Since we seek an 
upper limit on the size of the emitting region, we use the speed of light as the propagation speed. 
With this assumption, the At = 100 s variability (rise and decay) observed during flare J21_2 
corresponds to an upper limit on the radial size of the X-ray emission region of < 3.0 x 10 12 cm 

(2.5 Schwarzschild radii, Rs), based on propagation at the speed of light from the center to the 
edge of a spherical region, disregarding the light travel time across the region. This corresponds 
to an upper limit on the emitting volume V < (4/3)7rc 3 Af 3 and, in turn, on the magnetic energy, 
Ub < B 2 V/8it, where B is a function of the radial distance from the SMBH, R. 

The ion internal energy, which could amount to ~10 times the magnetic energy (e.g. Hirose et al. 
2004), is not taken into account here as it is not clear how it would be converted into photons on 
such a short timescale. In any case, including the ion energy would require us to use the plasma 
sound speed as the relevant propagation speed (significantly smaller than c), which would compen- 


- 16 


sate the effect of increased energy density by making the emitting volume much smaller and would 
lead to even tighter constraints on the flaring region location. 

We estimate the strength and radial profile for B based on theoretical considerations and previ- 
ous accretion disk simulation results. In the jet, a R" 1 dependence is predicted (Blandford & Konigl 
1979). For a disk, the magnetic field strength is estimated from equipartition of magnetic and 
kinetic energy. The theoretical dependence for a radiatively inefficient accretion flow (RIAF) 
is B oc R~ 105 and Bondi flows and Advection Dominated Accretion Flows give B oc R ~ 125 
(Yuan, Quataert & Narayan 2003; Narayan et al. 1998), while simulations typically give between 
R ~ and R ~ 1-2 (Dexter et al. 2010). Thus, we conclude that the disk profile is very close to 
the jet profile, and we only show one curve, with i? -1 , in Figure 8. Values in the inner part 
of the disk in the range of 10-100 G have been quoted in previous work (Falcke &; Markoff 2000; 

Yuan, Quataert & Narayan 2003; Dodds-Eden et al. 2010). Thus, we set the value to be Br s = 

100 G at 1-Rs) and the exact profile we assume is shown in Figure 8b. We view this as being a 
conservative assumption in the sense that a lower magnetic field would constrain the flare to occur 
even closer to the SMBH. 

Figure 8a shows the main result, where we plot the photon production efficiency, r /, which is 
the ratio of the fluence of the flare to the energy available at a given radial distance from the SMBH. 

As indicated by the fast rise and decay observed in Flare J21_2, we assume that the emitting region 
keeps a constant size during the flare, and use the full duration of the flare (1,898 s) and its average 
luminosity of 2.1 x 10 35 ergs -1 for this calculation. Thus, based on energetics, there is only enough 
energy to power the flare if it originates within < 10 i?s (Br s / 100 G) t ? 1 / 2 of the SMBH. 

Here, we discuss how beaming could affect this conclusion. If the flare emission is beamed 
towards us, the rest frame variability timescale would be larger than the observed one, and the 
actual luminosity would be lower than we derived. Beaming becomes significant when the bulk 
motion of the emitting particles has a relativistic speed in a direction within a few tens of degrees 
from our line of sight. However, we argue here that this case is unlikely, and physically reasonable 
levels of beaming would have only a mild effect. 

The spin axis of the SMBH is largely unconstrained, but recent studies suggest that the spin 
axis points >35° away from our line of sight (Dexter et al. 2010; Broderick et al. 2011; Shcherbakov, Penna &; McKi 
2012). In jet models, the X-ray emission comes predominately from the base of the jet, the so-called 
nozzle. The bulk speed in the nozzle is assumed to be the sound speed, ~ 0.4c (Falcke &; Markoff 
2000). The Lorentz factor 7 of the bulk motion in the nozzle is thus ~ 1.1, which yields a Doppler 
factor D = 1/(7 (1 — {3 cos 9)) of ~ 1.4 for the extreme case of 9 = 35°, very mildly affecting 
the timescales, and decreasing by a factor of 2 ( D 2 ) the emitted luminosity (see consequences on 
Figure 8a). 

Alternatively, one could think that a hot spot orbiting the SMBH could produce beamed 
emission if its orbital path were tangential to our line of sight. There is indeed evidence indicating 
that we may be seeing the disk nearly edge-on (Hamaus et al. 2009). However, typical orbital 
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velocities are in the 0.3 - 0.4 c range, which does not produce strong beaming. Moreover, this 
configuration would produce a modulation at the orbital period (between ~ 1860 s and ~ 360 s 
at the innermost stable circular orbit for a SMBH that is non-rotating and maximally rotating, 
respectively), which has not previously been seen in X-ray observations. 



2 4 6 8 10 12 14 16 

R/R s 


Fig. 8. — Constraining the location of the flare based on the energy available in the Sgr A* disk 
or jet. Panel a shows the required efficiency for converting the available magnetic energy into the 
photons we see in flare J21_2. Panel b shows our estimate (see text) for how the magnetic field 
strength ( B ) depends on the radial distance from the SMBH ( R ) in units of Schwarzschild radii 
( Rs = 2 GM/c 2 , where G is the gravitational constant, M is the mass of the SMBH, and c is the 
speed of light) . 


9. Summary 

Four flares were detected during this observation campaign, two of medium amplitude, J21_2 
and 017, and two weaker ones, J20 and J21_l. For the first time, we clearly see that the power-law 
spectrum of Sgr A* X-ray flares extends to high energy, with no evidence for a cut off; flare 017 is 
detected up to 79 keV and flare J21_2 up to 60 keV. 

The SED fit slightly favors the SB model for flare J21_2 but is inconclusive for flare 017 due 
to limited statistics at high energy. However, we confirm previous reports that physical arguments 
favor the SB model over the two simplified IC models we tested. This implies that efficient particle 
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acceleration takes place continuously over the duration of the flares, as the synchrotron cooling 
time is of the order of 1 s. 

Flare J21_2 was detected over a long enough time to search for spectral evolution, but we find 
none. We detect a variation of photon index between two individual flares at the 95% confidence 
level. Flare 017 (r = 2.04+q2q), which is also stronger, is harder than flare J21_l (r = 2.84 %q^). 
This is opposite to the finding of Degenaar et al. (2013) who reported on one bright flare that was 
softer than the average of five weak flares. Keeping in mind that the present result and that of 
Degenaar et al. (2013) have low significance, this indicates that the photon indices of flares can take 
a range of values, even for a given flux, which is suggestive of magnetic reconnection as acceleration 
mechanism. 

Variability with a timescale of 100 s was observed during flare J21_2, both in rise (factor of 
at least 3.8 ± 1.1) and decay (factor of 3.4 ± 1.8). This suggests that the flaring region keeps the 
same size throughout the flare, and it allows us to put an upper limit of 2.5 Rs on its radial size. 
Assuming that flares are powered by magnetic energy, this flaring region size places a constraint 
on the location of the region: energetics show that there is only enough energy to power flare J21_2 
within 10 Rs of the event horizon. 
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A. Bayesian block search sensitivity 

In order to characterize the sensitivity of the peak search using the Bayesian block method 
(section 3), a set of synthetic event lists featuring a 1,000 s long rate increase were produced. The 
baseline rate was chosen to be 0.6 cts/s and the peak rate was gradually increased from 0.62 cts/s 
to 0.8 cts/s by steps of 0.02 cts/s. Fifteen event lists were generated per peak count rate value, 
making a total of 150 generated sets. The Bayesian block search was then applied to each event list, 
and the success of the search was evaluated by looking at the presence and location of the change 
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points; Two change points, each within 1,000 s of the actual time of rise and decay is considered a 
successful search, even if more change points are detected. The actual peak significance is calculated 
using the knowledge of the actual rise and decay times, and 3,000 s before and 3,000 s after. It 
varies from set to set as a random generator is used to generate the intervals separating events, 
following a Poisson distribution. Finally, the trials are binned by actual peak significance, and for 
each bin, the ratio of the number of success to the number of trials is plotted. As we use a geometric 
prior for the Bayesian block search (Scargle et al. 2013), the sensitivity depends on the number of 
events. Figure 9 shows the results for 52,109 events and 28,000 events, which correspond to the 
July and October light curves, respectively. 

This study shows that our settings ( fpr = 0.01, and geometric prior described in section 3) 
yield a detection chance > 50% for a peak significance between 4.5 a and 5 a and of 100% for a 
peak significance >5 a. Performing these simulations with different peak duration leads to similar 
results, only the significance of the peak matters. So it appears that flare J20 detected during the 
July observation with a significance of 4.8 a is at the limit of detection. This corresponds to a 
strength of ~ 14 times the quiescent flux, for a flare of ~ 920 s. 
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Fig. 9. — Calibration of the Bayesian block peak search using synthetic event lists made of 52109 
events (red solid line), and 28000 (green dash line). This plot shows the probability of detection of 
a rate increase lasting 1000s within a constant rate of 0.6 cts/s. 


B. Lightcurve of flare J21_2 

Figure 10 shows again the 3-79 keV light curve of flare J21_2 but without combining the two 
focal plane modules. The bin size is set to 200 s in order to yield comparable statistics with Figures 
2 and 3. One can see that the sharp rise and sharp decay are present in both modules (see section 

4 ). 
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Fig. 10.-- Light curve of the J21_2 flare (on 2012 July 21) with 200 s time resolution and 1-u 
upper and lower count rate limit. FPM A (black crosses) and FPM B (blue plus signs) are not 
combined in this plot. The data is extracted from a 50" radius extraction region centered on the 
radio position of Sgr A* and is not background subtracted, as in Figures 2 and 3. The green vertical 
dash lines define the sub-flare (J21_2a), and the rise (r), and peak (p), and decay phases of flare 
21_2b, as in Figures 3 and 5. 


C. Flares flux tables 

The energy density for flares J21_2 and 017 are given in Tables 4 and 5, respectively. We used 
the following method to calculate these SEDs. We start from the unfolded photon flux density (in 
photons cm -2 s” 1 keV” 1 ), as returned by ISIS (Houck & Denicola 2000b), and apply the absorption 
correction. Then, the energy density points are obtained by multiplying by 1.602 x 10” 9 E gco , 
where E geo is the geometric mean of each energy bin, defined as E geo = \[E\E^. We use ISIS to 
unfold the spectrum i.e., to account for the detector response matrix and the effective area of the 
optics. For the absorption correction, we used the column density of 1.66lg' 6 l j ) x 10 23 cm” 2 that 
we determined with the joint fit of the four flares. The 90% errors are divided by a factor of 1.7 
to obtain the 1-a range Ah = 1.30-2.07x 10 23 cm” 2 . For all three values (1.30 x 10 23 , 1.66 x 10 23 , 
and 2.07 x 10 23 cm” 2 ), we calculated the ratio of unabsorbed to absorbed flux, and we multiplied 
the absorbed flux values by these ratios to make the absorption correction and to determine the 
errors on the absorption correction. 
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Table 4. Spectral Energy Distribution for flare J21_2. 


E (keV) 

AE/2 (keV) 

E Fe (erg cm 2 s x ) 

Error (erg cm 2 s x ) 

Significance (a) 

3.68 

0.68 

1.50314e-ll 

4.04299e-12 

7.78 

4.92 

0.52 

1.04005e-ll 

1.73769e-12 

7.36 

5.90 

0.42 

1.16285e-ll 

1.68437e-12 

7.48 

7.02 

0.66 

1.04869e-ll 

1.52881e-12 

7.09 

8.72 

1.00 

8.27270e-12 

1.17714e-12 

7.29 

10.94 

1.18 

9.54882e-12 

1.31518e-12 

7.29 

14.36 

2.20 

1.02458e-ll 

1.42814e-12 

7.45 

19.98 

3.38 

7.90110e-12 

1.74097e-12 

4.55 

41.70 

18.30 

7.23571e-12 

2.25508e-12 

3.02 


Note. — The first column gives the midpoint energy of the bin, the second, the half-width of 
the bin, the third, the specific flux multiplied by E geo , and the fourth, the l-cr error (including 
the propagation of the uncertainty on the column density). 


Table 5. Spectral Energy Distribution for flare 017. 


E (keV) 

AE/2 (keV) 

E Fe (erg cm 2 s 1 ) 

Error (erg cm 2 s 1 ) 

Significance (a) 

3.82 

0.82 

2.05094e-ll 

5.12969e-12 

7.70 

5.14 

0.46 

1.81119e-ll 

2.86282e-12 

7.55 

6.38 

0.74 

1.48637e-ll 

2.11914e-12 

7.44 

7.88 

0.72 

1.55779e-ll 

2.26569e-12 

7.40 

9.50 

0.86 

1.65479e-ll 

2.25693e-12 

7.55 

11.60 

1.20 

1.80809e-ll 

2.47585e-12 

7.38 

14.88 

2.04 

1.73249e-ll 

2.64067e-12 

6.59 

19.28 

2.32 

2.06093e-ll 

3.82042e-12 

5.42 

50.32 

28.68 

9.99168e-12 

3.36686e-12 

2.98 


Note. — The columns description is identical to that in Table 4 
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